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Abstract. Given u > i and <5 > arbitrary, we prove the existence of energy 
solutions of 

(0.1) d tt u-Au~u 5 =0 

in R 3+1 that blow up exactly atr = i = 0ast-»0- . These solutions are 
radial and of the form u = X(t) 2 W(X(t)r) + rj(r, t) inside the cone r < t, where 
X(t) = t" x ~" , W(r) = (1 + r 2 /3)~i is the stationary solution of JOT}, and T} 
is a radiation term with 



/ 

J\r 



(|Vf?(ai, t)\ 2 + \r) t (x, t)\ 2 + \ v (x, t)\ 6 ) dx^O, t->0 

r<t] 

Outside of the light-cone there is the energy bound 

/ (|Vu(x,t)| 2 + \u t (x,t)\ 2 + |u(a:,i)| 6 ) dx < 5 
J[r>t] 

for all small t > 0. The regularity of u increases with v. As in our accompa- 
nying paper on wave-maps [101 . the argument is based on a renormalization 
method for the 'soliton profile' W(r). 

1. Introduction 

Since the seminal paper of Jorgens [5] much work has been devoted to the study 
of well-posedness of the nonlinear wave equation 

d tt u - Au + f(u) = 

in R 3 ^ 1 and suitable nonlinearities f(u). Jorgens showed that for H 1 (R 3 ) subcriti- 
cal defocusing nonlinearities f(u) = |u| p_1 u withp < 5 smooth data lead to smooth 
solutions for all times. The critical defocusing casep = 5 was resolved by Struwe [H] 
for radial data and Grillakis [S] for general data. These authors proved global well- 
posedness and scattering results for energy solutions, see Shatah-Struwe [17] and 
Sogge [15] . No corresponding results are known for the supercritical case p > 5. 

In this paper we address the solvability of the nonlinear wave equation in M 3+1 
with a focusing nonlinearity f(u) = — |u| p_1 m. In this case blow-up may occur. 
Indeed, it was shown by Levine [llj via a convexity argument that data in (H 1 n 
L 2 )xL 2 with negative energy lead to finite-time blow-up, see also Strauss [50] . Local 
well-posedness in the optimal regularity class was considered by several authors, see 
Sogge [18] for a detailed exposition of this work. Most relevant for us is the case 
p = 5 where the equation is locally well-posed in the energy space H 1 x L 2 (R 3 ). 
Moreover, if the solution cannot be continued beyond some finite time T* as an 
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energy solution, then necessarily the Strichartz norm ||u||i8([ 0! T,)xR 3 ) — 00 (with 
similar results in all dimensions). 

The question of the blow-up rate was addressed by Merle-Zaag in the conformal 
range p < 3, see jT3]— [TS] (their results extend to all dimensions). They showed 
that if solutions to the Cauchy problem, with 1 < p < 3, 

d tt u - Au ~ = 0, (it , ui) G Hl oc x Lf oc 

blow up in finite time T*, then the following holds: for any a£l 3 the self-similar 
change of variables 

u(x, t) = (T* - ty^w a {y, s), y = ^ — - s = - log(T* - t) 

leads to functions w a satisfying 

sup \\w a (s)\\ H i(B) + \\d s w a (s)\\ L 2 {B) < K 

s>- logT,+l, aeM 3 



where B is the unit ball and a constant K that only depends on p, T* and the norm 
of the initial data in Hl oc x Lf oc . 

For the energy critical case p — 5, i.e., 



(1.1) d tt u-Au-u b = 

there has been some recent activity, see [9] , [8] , [7] , which we now describe in more 
detail. The Talenti-Aubin solutions 

W(r) = (1 + r 2 /3)-^ 

are extremizers of the Sobolev imbedding if 1 (R 3 ) =— ► L 6 (M 3 ) and satisfy the Euler- 
Lagrange equation —AW — W 5 = 0. In [5] the first two authors showed that there 
exists a small co-dimension one manifold M. around W in a suitable topology so 
that data on this manifold exhibit global existence and an asymptotic behavior of 
bulk term plus radiation. The radiation term is also shown to scatter like a free 
energy wave. It is conjectured, see [2], that this manifold has the property that 
it separates a region of scattering from one of blow-up. As a first result in this 
direction, Karageorgis-Strauss [7] showed that above the tangent space of M. at W 
finite time blow-up occurs, albeit for the equation 

d tt u-Au- \u\ 5 = 

Note that the result of [9j equally well applies to the nonlinearity |u| 5 (in fact, the 
solutions constructed in [9] are positive so that there is no distinction between u b 
and I it 1 5 from the point of view of that paper). 

Kenig-Merle [8] studied the behavior of solutions with data {uq,Ui) G H 1 x 
L 2 (M 3 ) of energy S(uq,Ui) < £(W, 0) where the conserved energy is 



S(u,u t )= ( [i(^ + |Vu| 2 )-^' 
J R 3 L2 6 . 



d.r 



They found that in this regime there is a dichotomy between blow-up and global ex- 
istence/scattering depending on whether ||Vuo||2 > ||VW||2 or ||V«o||2 < HVWlj 2- 
Note that 

W(x,X) := X^W(Xx) 
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is a stationary solution of (|1.1|) for all A > 0. Moreover, the energy is constant in 
A (reflecting the energy criticality of the equation). Linearizing the wave equation 
around W leads to the linearized operator 

H = -A - 5W 4 

The wave evolution of H has two types of instabilities: an exponential instability 
arising from the negative spectrum of H (which has a unique negative eigenvalue) 
as well as a "bound state" at zero energy: H(d\W\\=i) — where d\W decays 
like r _1 and thus does not belong to L 2 (K. 3 ) — this is what one refers to as a zero 
energy resonance. In this paper we construct blow-up solutions by 'projecting out' 
the exponentially growing mode of the linearized equation. 
More precisely, we seek radial, real- valued, blow-up solutions 

u(x,t) = A(f)*W(A(t)x) +r)(x,t) 

of where A(t) — > oo as t — > and with the local energy inside the light-cone 
\x\ < t of T](x,t) going to zero as t — > 0. The local energy relative to the origin is 
defined as 

£ loc (v)= I (v? + \Vv\ 2 + \v\ 6 )dx 

J[\x\<t] 

The following theorem is the main result of this paper. The blow-up occurs at time 
t = when solving backwards in time. 

Theorem 1.1. Let v > ^ and S > 0. Then there exists an energy solution u of 
(jl.ip which blows up precisely at r = t = and which has the following property: 
in the cone \x\ = r < t and for small times t the solution has the form, with 
\{t)=t- x - v , 

t) = \i(t)W(\(t)r)+ V (x,t) 
where £\ oc (r](-,t)) — > as t — > and outside the cone u(x,t) satisfies 

[ [\Vu(x, t)\ 2 + \u t (x, t)\ 2 + \u(x, t)\ 6 ] dx < 5 

J[\x\>t] 

for all sufficiently small t > 0. In particular, the energy of these blow-up solutions 
can be chosen arbitrarily close to £{W, 0), i.e., the energy of the stationary solution. 

The restriction v > 1/2 arises only due to technical reasons, and we hope to 
eliminate it in subsequent work. If v > 1, then the solutions from Theorem 11.11 
belong to L°°(R 3 ) for alH > and blow up at the rate 

IM,f)||oo*H 1+ "V2 

as t — > 0. The proof is based on a renormalization procedure analogous to the 
one that the authors used for the construction of blow-up solutions for wave maps 
in [10] . For our purposes this refers to the fact that we do not simply perturb around 
X^(t)W(X(t)r) to obtain the linearized equation for -q, but rather first modify the 
blow-up profile and then perturb around this " rcnormalizcd" profile. More precisely, 
fix a large integer N. Then there exists a function u e satisfying 

(1.2) u e e C*¥-({to > t > 0, |x| < £}), £i oc (u e ){t) < (tX(t))- 1 as t ^ 
so that the radiation term rj above has the form 

rj(x,t) = u e (r,t) +e(r,t), 0<r<t 
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where e decays at t = 0. In fact, r\ can be extended globally with the property that 
£e^-(l 3 ), s t € ^-^-(M 3 ), £ioc(e)(*) <** as ^0 

with spatial norms that are uniformly controlled as t — » 0. 

As this paper was written concurrently with our wave-map paper |10] it is only 
natural that there would be some similarities between this paper and [10] . In fact, 
a secondary goal here is to show that the method used in both papers is flexible 
and applies to quite distinct scenarios. The main differences between this paper 
and [10] are as follows: 

• The blow-up profile is not constant in L°° but rather grows at rate A 3 . The 
renormalization procedure thus needs to be adapted to this case. 

• In contrast to [10] , the linearized operator exhibits negative spectrum. This 
produces exponential instability of the linearized wave flow. 

• The linearized operator no longer exhibits a strongly singular potential in 
the sense of [4]. Thus, a (Dirichlet) boundary condition is needed at R = 0. 

We feel that the most important difference listed here is the exponential instabil- 
ity. In fact, as in the asymptotic stability paper [9], we need to 'project out' this 
exponential growth. Our blow-up rates are therefore expected to be non-generic. 

2. The renormalization step 

In this section we show how to construct an arbitrarily good approximate radial 
solution to the wave equation (|1.1[) as a perturbation of a time-dependent ground 
state profile 

u = X^W(R), W(R) = {1 + R 2 /3rK R = rX(t) 
with the polynomial timescale 

A(t) = r 1 -", v > 0. 

Theorem 2.1. Let k G N. There exists an approximate solution U2k-\ for 
of the form 

u 2k ^(r, t) = As ( t ) [W(R) + ^i? 2 (l + i?T 1/2 + O ( ^ 2(1 ( +jr r 
so that the corresponding error has size 

e2k~ 



1 yt 2 {txf k j 



Here the O(-) terms are uniform in < r < t and < t < to where to is a fixed 
small constant. 

Remark 2.2. The u e in Q is 

The analysis below shows that it has the stated regularity up to the light-cone. 
Moreover, one checks that 

A^ 
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which is the claimed decay rate for the local kinetic energy of u e . The local potential 
energy of u e decays like (iA)~ 3 . 

Proof. We iteratively construct a sequence u k of better approximate solutions by 
adding corrections v k , 

Uk = v k + Ufc-i 

The error at step k is 

2 

e k = {-d 2 + d 2 + -d r )u k + u% 
r 

If u were an exact solution, then the difference 

e = u — tifc-i 

would solve the equation 
2 

(2.1) (-9 4 2 + ^ + -a r )e + 5^_ 1 £ + 10u|_ 1 e 2 + 10^_ 1 e 3 + 5ii fe _ie 4 +e 5 + e fc _i = 

r ..... 

In a first approximation we linearize this equation around e = and substitute 
Uk-i by uq. Then we obtain the linear approximate equation 

(2.2) (^d 2 + d 2 + ^d r + 5u{\ e + e fc _! « 

For r<fwe expect the time derivative to play a lesser role so we neglect it and 
we are left with an elliptic equation with respect to the variable r, 

(2.3) ^a 2 + ^ I . + 5u^e + e fe _i«0, r<£ 

For r w t we can approximate Up by zero and rewrite (|2.2p in the form 

(2.4) ^-d 2 + a 2 + ^e + e fc _ 1 «0 

Here the time and spatial derivatives have the same strength. However, we can 
identify another principal variable, namely a — r/t and think of e as a function of 
(i, a). Later, we reduce the above equation to a Sturm-Liouville problem in a which 
becomes singular at a = 1. 

The above heuristics lead us to a two step iterative construction of the w^'s. The 
two steps successively improve the error in the two regions r <gC t, respectively r«t 
To be precise, we define v k by 

(2.5) (d 2 r + \ + 5u£j v 2k+ i + e° 2k = 
respectively 

(2.6) (-d? + d 2 + ^d r ^V2k+e° 2k _ 1 =0 

both equations having zero Cauchy dat£0 at r = 0. Here at each stage the error 
term e k is split into a principal part and a higher order term (to be made precise 
below) , 

e k = el + el 



lr The coefficients are singular at r = 0, therefore this has to be given a suitable interpretation 
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The successive errors are then computed as 

e 2k = e\k-\ + N 2k{v2k), e 2 fe+i = e\ k - d 2 v 2k+x + N 2k+ i(v 2k+ i) 

where 

(2.7) N 2k+1 (v) = 5(u 4 2k - -4) v + 10u 3 2k v 2 + I0u 2 2k v 3 + hu 2k v 4 + v 5 
respectively 

(2.8) N 2k (v) = f>u\ k _ r v + lOul^v 2 + 10u^_ 1 « 3 + hu 2k ^v 4 + v 5 

To formalize this scheme we need to introduce suitable function spaces in the cone 

C = {{r,t) : 0<r<t, < t < t Q } 

for the successive corrections and errors. We first consider the a dependence. For 
the corrections Vk we set 

v — 1 1 

and use 

Definition 2.3. For i G N we let = if is is irrational, respectively j(i) = i if 
v is rational. 

a) For any positive integer k, we define Q to be the algebra of continuous func- 
tions q : [0, 1] — ► K with the following properties: 

(i) q is analytic in [0, 1) with an even expansion at 0. 

(ii) Near a = 1 we have an absolutely convergent expansion of the form 

oo 

q(a) = q (a) + £(1 - af^~ 2 ^ £ ««(o)(Iog(l - a)y 

i=l ]=0 

with analytic coefficients q , qij. 

b) Q m is the algebra which is defined similarly, with the additional requirement 
that 

3ij(l) = if i = 4fc + l > 4m+l. 

We remark that the exponents of 1 — a in the above series are all positive because 
of Po > —\- For the errors we introduce 

Definition 2.4. Q' is the space of continuous functions q : [0, 1) — > M. with the 
following properties: 

(i) q is analytic in [0, 1) with an even expansion at 

(ii) Near a = 1 we have a convergent expansion of the form 

oo 

q(a) = q (a) + ]T(1 - af^- 2 ^ 1 ]T gy (a)(log(l - a)Y 

i=l 3=0 

with analytic coefficients qo, q^ . 

b) Q' m is the space which is defined similarly, with the additional requirement 
that 

q^ (1) = if i = 4k + 1 > 4m + 1. 
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By construction, Qk C Q' k . The families Q! and Q' k are obtained by applying 
a~ 1 d a to the algebras Q and Qk, respectively. 

We remark that the number of logarithms in these definitions in the case when 
v is rational is far from optimal, but we have chosen this form since it simplifies 
the presentation. Next we define the class of functions of R: 

Definition 2.5. S m (R k ) is the class of analytic functions v : [0,oo) — ► K with the 
following properties: 

(i) v vanishes of order m and R~ m v has an even Taylor expansion at R = 0. 

(ii) v has a convergent expansion near R — oo, 



k-2i 



i=0 



The importance of even expansions in R lies with the fact that only those corre- 
spond to smooth functions in R 3 . For the same reason, we will work with even m. 
We also introduce another auxiliary variable, 



^ 6 =w 

Since we seek solutions inside the cone we can restrict b to a small interval [0, bo]- 
We combine these three components in order to obtain the full function class which 
we need: 



Definition 2.6. a) S m (R k , Q n ) is the class of analytic functions v : [0, oo) x [0, 1] x 

[0, bo] -> K so that 

(i) v is analytic as a function of R,b, 

v : [0, oo) x [0, bo] -> Qn 

(ii) v vanishes of order m and R~ m v has an even Taylor expansion at R = 0. 
(Hi) v has a convergent expansion at R — oo, 



v(R,;b) = Y,Ci(;b)R k ~ 2i 

i=0 

where the coefficients Ci : [0, bo] — > Q m are analytic with respect to b. 

b) lS m (R k , Q n ) is the class of analytic functions w on the cone Co which can be 
represented as 

w(r,t)=v{R,a,b), veS m (R k ,Q n ) 

We note that the representation of functions on the cone as in part (b) is in 
general not unique since R, a, b are dependent variables. Later we shall exploit this 
fact and switch from one representation to another as needed. We shall prove by 
induction that the successive corrections Vk and the corresponding error terms 
can be chosen with the following properties: For each k > 1, 
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(2.10) V 2k -1 G T^fe K 2 (^, Qk-l) 



(2.H) t 2 e 2k ^ 1 e 77 ^lS Q (R,Q' k _ 1 ) 



(txy 
(txy 



(2-12) ^G^y^IS 2 ^ 3 ^) 

(2.13) t 2 e 2k G j^k-[lS (R-\Q k ) + blS (R,Q' k )] 

with (|2.13[) also valid for k = 0. We remark that the order of vanishing at R = 
can be successively improved with k, but this does not appear to be important. 

Step 0: The analysis at k = 

With our notations, one checks that 

(2.14) t 2 e = -t 2 dtt[X^W{X(t)r)] G X^IS°(R- 1 ) 

as claimed. Now assume we know the above relations hold up to k — 1 with k > 1, 
and we show how to construct v 2k -i, respectively v 2kl so that they hold for the 
index k. 

Step 1: Begin with e 2k - 2 satisfying (|2.13[) or (|2.14[) and choose v 2k _\ so that 
(f2~T0l) holds. 

If fc = 1, then define e[J := eo- If k > 1, we use (|2.13|) to write 

e2fc-2 = e 2fe _2 + e 2fc-2 

where 

* e K-2 G 7^Z2 IS°(iT \ Q fc _x), i 2 e 2fe _ 2 G ^fe IS°(iZ, Q' fe _ x ) 

In the first term we can set b — and eliminate the b dependence, as all the b 
dependent part can be included in the second term. 

We note that the term e 2fc _ 2 can be included in e 2 fc-i, cf. (|2. 1 1[) . We define 
v 2k -i as in (|2.5[) neglecting the a dependence of e 2fc _ 2 . In other words, a is treated 
as a parameter. Changing variables to R in (|2.5p we need to solve the equation 

(tX) 2 Lv 2k ^ = t 2 e° 2k _ 2 G fc _ 2 IS ^" 1 , Q fc _x) 
where the operator L is given by 

L = -d 2 R - ld R - 5W 4 (R) 
Then (|2. 10|) is a consequence of the following ODE lemma. 
Lemma 2.7. The solution v to the equation 

Lv = / G S°{R- 1 ), v(0) = v'(0) = 

has the regularity 

v G S 2 (R) 
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Proof. Since / is analytic at with a constant leading term, one can easily write 
down an even Taylor series for v at with a quadratic leading term. 

It remains to determine the asymptotic behavior of v at infinity. For this it is 
convenient to remove the first order derivative in L (to achieve constancy of the 
Wronskian). Thus, we seek a solution of 

LR v = Rf, L = d% + 5W 4 = d 2 R 



(l + i? 2 /3) 2 

We use this fundamental system of solutions for L: 

4>{R) = R(l - i? 2 /3)(l + i? 2 /3)"3 

9(R) = (1 + i? 2 /3)^(l - 2R 2 + i? 4 /9) 

Clearly, Ld\W = and we set <f> = Rd\W\\ = \ up to a constant. The function 
9 is then determined from the Wronskian constancy condition W(6, cf>) = 1. This 
allows us to obtain an integral representation for v using the variation of parameters 
formula, which gives 

[■R rR. 

v = —R~ 1 9(R) / cj)(R')R' f(R!) dR' + RT^R) / 9(R')R'f(R') dR' 
Jo Jo 

The right-hand side grows like R, as claimed. □ 



(2.15) Vl= v^R), VeS 2 (R) 



As a special case of the above computation we note the representation for v\ , 

(tXf 

This justifies the choice of the second term in the expansion for it2fc-i in Theo- 
rem [2TTJ 

Step 2: Show that if V2k-i *s chosen as above then (|2.1ip holds. 

Thinking of i^fc-i as a function of t, R and a we can write e2k-i in the form 

e 2 fe-i = N 2 k-i(v 2 k-i) + E t v 2 k-\ + E a v 2k -i 

Here N 2 k-i(v2k-i) accounts for the contribution from the nonlinearity and is given 
by (|2.7[) . E t V2k-\ contains the terms in 

(2.16) -d tt v 2 k-i(t,R,a) 

where no derivative applies to the variable a, while E a V2k-i contains those terms 



(d tt - rr - ^dr^j U 2 fc-1 (*, R, a) 



r 

where at least one derivative applies to the variable a (recall that in Step 1 the 
parameter a was frozen). We begin with the terms in Nzk—i- We first note that, 
by summing the Vj over 1 < j < 2k — 2, 

(2.17) u2 k -2-u £ -^lS 2 (R,Qk-i) 

(txy 

The first term in iVafc— i(u2fe— l) contributes 

t 2 ( u 2k-2 ~ U o) V 2k-l = < 2 [(u 2 fe-2 - Mo) 4 + 4(u 2 fe-2 - «o) 3 Wo 

(2.18) +6(u 2 fc- 2 - u ) 2 ul + 4(u 2/£ - 2 - u )ul}v 2 k-i 
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Using (|2.17|) we compute 



A^ 



C a 6 IS 2 (iT 2 , Q fe -i)77^ IS Z (-R, Q fe _i) 
A* 



c (^F Is2(irl ' Sfc - 

as well as 



t 2 {u 2k -2 - uo)u 3 v 2 k-i G t 2 ^ IS 2 (i?, Q fc _ 1 )Ai5°(i?- 3 )— ^IS 2 ^, Q fe -i) 



The other two terms in (|2.18p are similar. Next, compute 



and 



A 



with similar statements for W2fe-2 u 2fe-i an d M 2fc-2f2fe-i- Summing up we obtain 

(|^ IS2(rl ' St - l)C (| 

This concludes the analysis of A^fc-i^fc-i)- We continue with the terms in 
E t V2k-i, where we can neglect the a dependence. Therefore, it suffices to com- 
pute 

^((^ IS2(i?) ) C (^ IS2(i?) 
Finally, we consider the terms in E a V2k-i- With 

(U) 2fe 

we have 

A 5 \ N A* 

(OP J «-a(i?,«) + OT 

-2Ra~ 1 w Ra - 2aT 1 w a (i?,a) + (a 2 - l)w aa (R, a) + 2aw a (R,a) 



V2k-i(r,t) = , w(R,a), w G 5 (i?, Qk-i) 



t 2 E a v 2k -i = -2tdt ( y-^2T aw a (i?, a) + — ^ [2(i/ + l)aJ?«; afl ( J R, a) 
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Since Qk-i are even in a we conclude that 

(1 - a 2 )d aa , ad a , a~ 1 d a : Qk-i -> Q^-i 

and therefore 

\| 

i^ 2fc _ lG _^IS 2 (i*, Ci_x) 

This concludes the proof of (|2.11|) . We remark that for the special case of k = 1, 
i.e., with «i as in (|2.15p . these arguments yield 

(2.19) t 2 ei G j^IS 2 (R) 



Step 3: Define v 2 k so that (|2.12[) holds. 

We begin the analysis with e2k-i replaced by its main asymptotic component e° 
around R = oo. This has the form 



which we rewrite as 



(2-20) t 2 4 k _, = 7^p9(a), q e QLi 



A^ 

t 2 e»-i = (a)2fc -i "g( a ) 
We remark that (|2.19[) implies that 

* 2 e~?(a) = ^a. 

Consider the equation (|2.6p with eS^^ on the right-hand side, 

2, 
r 

We look for a solution fi2fc which has the form 

A^ 



t 2 ( -d 2 + d 2 + -8 r i 2k = -#e\ h _ x 



V2k 



Wr~' W2k{a) 



Thus, 

t 2 [-d 2 + d 2 + h^j U^W 2k (a)^ = —^aq(a) 
Conjugating out the power of t we get 

t 2 ( - (ft + ( 2fc -l>-^ 2 + 9 2 + Hft) ^ (o) = ag(a ) 

which we rewrite as an equation in the a variable, 

(2.21) L (2fc-!)»-| W M = aq(a), q G 
where the one-parameter family of operators Lp is defined by 

(2.22) Lp = (1 - a 2 )d aa + 2(a- 1 + a/3 - a)9 a - (3 2 + f3 

We claim that solving this equation with zero Cauchy data at a = yields a solution 
which satisfies 

(2.23) W 2 k{a) = a 3 qi (a), q x G Q fe 
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This gives 



v * k = (t\)2k-i a 91 (°) = ( tA )2fc+2 J? 91 (a) 



which is not entirely suitable as the next correction since it has an odd expansion 
at R — instead of an even one. To remedy this we simply set 

- R3 (irW 9l(a) e (^ 54(i?3 ' Sfe) 

Clearly, this will conclude Step 3. To prove the claim (|2.23p we need the following 

Lemma 2.8. Let f G aQ' k l , k > 1. Then there is a unique solution w £ a 3 Qk to 
the equation 

(2.24) L (2k _s )v _iw = f, w(0) = 0, d a w(0) = 
Proof. Denote 

(2.25) /? = (2fc - - 1 = (4fc - 3)/3 + 2(fc - 1) > ~. 
We write 

= a- 2 d a (a 2 d a ) - a 2 9 QQ + 2(/3 - l)a9 - /3 2 + f3 

To study the behavior of the solutions at we match the coefficients in Lp w = f 
with 

oo oo 
3=1 3=2 

yields the system, with j > 1, 

(2j + l)(2j + 2)^+1 + [-(2j - l)(2j - 2) + 2(/3 - l)(2j - 1) + (/? - /3 2 )] Wj = 

where we take i«i = 0. The coefficient of Wj+\ is always nonzero; this allows us 
to successively compute the coefficients uij. The convergence of the series for w 
follows from the convergence of the series for /. 

It remains to study the solution w near a — I. The behavior of Lp at 1 is well 
approximated by 

L} 3 = 2(1 - a)d aa + 2(3d a = 2(1 - a)" +1 o [(l - a)~ fi d a ] 

which annihilates the functions 1 and (1 — a)^ +1 . Therefore, we seek a fundamental 
system for Lpy — of the form 

oo oo 

(2.26) &(o) = l + ^/i,(l-a) f , 2 (a) = (1 - af +1 [l + ^ fi e (1 - a) 1 
This leads to the conditions, with no = (Xq = 1, 

t-i 

(2.27) 2nt +1 (t +l)(t-p) = hW -2(3 + 3) + (3 2 -P}+ 2 £ j H 

3=0 

(2.28) 2fi e+1 (e + (3 + 2)(£+l) = fi e [([3 + I + £)(£ - [3 + 4) + [3 2 - [3] 

i-\ 

+ 2^ + 1+^ 
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Clearly, (|2.28|) always has a solution whereas (|2.27p requires (3 £ Zj; in the latter 
case, the series in (|2.26p define entire functions. If, on the other hand, (3 £ Zq , 
then cf>i is modified to 

oo 

(2.29) <t>i{a) = 1 + ^^(1 -a) £ + Cl 2 (a)log(l -a) 

e=i 

with some unique choice of c\. 

Modulo a linear combination of 4>\ , 4>2 it suffices to find one solution to the 
inhomogeneous equation Lpw = f near a = 1. We begin with the case when / has 
the form 



/( ) = (i- a y>XA(i- 

with 7 > 0. Then we seek a solution w of the form 



a) 



w{a) 

k=0 



This leads to the system 
(2.30) 

2^o(7+l) (7-/3) =/o 

2w e+1 ( 1 + i + 2)( 7 -13 + 1+1) =f t + «j € [( 7 + 1 + 1)( 7 + 1 + 4 - 2/3) + /3 2 - (3\ 



2e£S(7+j + i: 



<7=0V i T J T 

which is solvable unless (3 — 7 is a nonnegative integer, in which case the solution 
w has the modified form 

00 

w(a) = (1 - a) 7+1 X>fc(l - + c0 2 (a) log(l - a) 

Indeed, assume that j3— 7 = p > with integer p. Then (|2.30[) is in general violated 
for 1+1 = p regardless of the choice of w p ; thus, c above has to be chosen so that 
the coefficients of (1 — a) 7+p = (1 — a)' 3 in Lpw — f match. This can be done since 

<h(a) log(l - a) = [(1 - af +1 + 0((1 - af+ 2 )} log(l - a) 

and 

L/s[0a(o) log(l - a)] = 2(/3 + 1)(1 - a)' 3 + 0((1 - a)' 3 "*' 1 ) 

with (3 + 1 > i and thus nonzero. Note that w p remains undetermined - this 
amounts to the freedom of adding a constant multiple of 02(a) to w(a). 
Similarly, if / has the form 

3 j 00 

/(a) = J2 (M 1 - a)) m fm(a) = ]T (log(l - a)) m (l - a)~< £ f km (l - a) k 

m=0 m=0 fe=0 

then we seek a solution u> of the form 

3 3 00 

w(a) = J2 ( lo g(! ~ a)) m Wm(a) = £ (log(l - a)) m (l - a) 7+1 £ «, fem (l - a) k 

m=0 m=0 fc=0 

Identifying the coefficients of the powers of log(l — a) we obtain the system 
L Wj = fj, Lpwj-x = fj-i - jQiwj, 
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respectively 

LpWj^k = fj_i - jQiWj-k+i + j(j - l)(9 2 w-,- fc+ 2, k > 2 

where 

l + o 2(1 + a) 2a(3 „ 1 + a 
1 — a a 1— a 1— a 

This system is solved iteratively as in the first case provided that (3 — 7 is not a 
nonnegative integer. Otherwise, the solution w has the modified form 

j 00 
w(a) = £ (log(l - a)) m [c m 2 (a) log(l - a) + (1 - a) 7+1 £ w km (l - a) k ] 

m=0 k=0 



Consider now / 6 Q'u-i- If v i s irrational then according to Definition 12.41 we 
can represent it in the form 

f(a) = f (a) + £(1 - aT^-^^Ma) 
i=l 

with fi analytic. Hence the exponent 7 above takes the values 

~i - 1 



7, = Ufa + 1) - 2 



- 1 



On the other hand, we have 



P = (2k - |)„ - 1 = (4k - 3)(A) + 1) - 2k + 1 = 7 4fe-3 

Then 7^-/3 can only be an integer if i = Ak — 3. However, in this case there is 
no logarithmic term due to the additional condition Q4fc_3(l) = which has the 
effect of replacing 74^-3 by 74/C-3 + 1. The conclusion of the lemma follows in the 
irrational case. 

On the other hand, if v is rational then / has the representation 

00 i 

f(a) = f (a) + £(1 - )*CA.+i)-*[*ii]-i £ / tf (o)(log(l - a))' 

i=l j=0 

with fij analytic. In order for — "fe to be a nonnegative integer we need to have 
I < Ak — 3. Again the condition 54^-3(1) = guarantees that there is no extra 
logarithmic term arrising from the i = Ak — 3 component of /. On the other hand, 
if I < Ak — 3 then we can contribute an extra logarithm for a total of I + 1 < Ak — 3 
logarithms to the i = Ak — 3 term of w. The conclusion of the lemma again follows. 

It is also worth discussing the special case k = 1. This will also serve to explain 
how the algebra Qk arises in the iteration. If k = 1, then (|2.21j) reduces to the 
equation, with the usual /3 = (v — l)/2, 

L/3o W 2 (a) = a 

due to t 2 fi(a) = \i(t\)~ 1 a. As discussed above, 
W 2 (a) = g (a) + 9l (a)(l - af« +1 if (3 & Z+ 

W 2 (a) = h (a) + hi(a)(l - a) r3o+1 + h 2 (a)(l - af 0+1 log(l - a) if O G 
Thus, we see that in all cases W 2 G Q\ for j = 0, 1 and a near 1. □ 
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Step 4: With v 2k as above show that e 2k is as claimed. 
Modifying (|2.20[) to insure an even expansion at R = we set 

2fc 1 (t\) 2k (R2 + l)i k 1 

Then we can write e2/c in the form 



k-l) + ^ (4k-! - (-9t + + 2 -d r ^j V 2k \ + f N 2k (v 2k ) 



t 2 e 2k = t 2 (e 2 fc_i - e 

where A^/c is defined by 

We begin with the first term in e 2k , which has the form 

t 2 (e 2fe _a - 4 k _ t ) G IS°(7r\ e' k _i) 

We claim that 

(2.31) IS°(i?- 1 ,Q' fc _ 1 )c IS°(i?- 1 ,Q fe _ 1 ) + ^ I IS 2 (i?,Q^ 1 ) 

For iu G IS^iZ -1 , Qfc_ x ) wc write 



tu = (1 — a 2 )w + R 2 w 

{txy 

Then 



as desired. 

Next, we consider the expression 

/ = * 2 (^ fc -i-(-d 2 + d 2 + V) « 2fe 

By construction this would vanish if the factor R(l + i? 2 ) _ 5 were dropped in both 
e 2fc-i an d w 2fc- Hence / contains only the components of the second term for which 
at least one derivative falls on the factor R(l + R 2 )~^: 

2„\ R „9„_„ R 



f = V 2k t Z ( -d 2 + 31 + -d r ) - ^ - 2t 2 d t V 2k d t 



(l+i? 2 )3 (l+i? 2 )2 



2t-d r v 2k d r - — T +v 2k — o r - 



(l + i? 2 )2 r (l + i?2) 5 

where t>2fc = v 2 kR{^ + R 2 )~^ ■ Then a direct computation gives 

as needed. 

Finally we consider the nonlinear terms in N 2k {v 2k ). Again the a, 6 dependence 
is uninteresting since Q k is an algebra. We start with the last term in (|2.8p . By 
construction, 

A3 _ . 
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Thus, 
Using that 



t 2 v 5 2k g IS 5 (i? 5 , Q k ) c {a)2k (u)8fc _ 2 (l + i? 2 ) 3 IS 1 ^" 1 , C*) 



(U) 8fc - 2 

we conclude that 



(1 + R 2 f _ ^ + + 3a 4 6 4fc-3 + a 6 6 4fe-4 



t 2 v* k €-^ IS 1 (R-\Q k ) 



which is an admissible contribution to (|2 . 13|) . Finally, we check the first term in 
dUg]) . From 

A 

(txy 

and the form of V2k, 

X^ „o,„ A 4 A* 



w 2 fe-i -wo£ 77777 IS 2 (i?, Qfe) 



^ II (/,') + -—IS^Q^J ' 2fc+2 IS 2 (R\Q k ) 



C ^A^iT 1 ) + X^a 2 IS°(R-\ Qk)) i IS 2 (i? 3 , G fc ) 



c (a) 2 is°(i?- 4 , g fc ) * 2 2 fc+2 IS"! A ) 



The other three terms in (|2.8p can be checked similarly. This concludes the proof 
of Theorem O □ 

3. The linearized problem 

We seek a radial solution of fTTTj) of the form it = u 2 fc-i + e with U2k-i as in 
Theorem 12. II Our ansatz leads to 



(3.1) d tt e -Ae- bX 2 (t)W A {X{t)x)e = N 2k -i(e) + e 2k -i 

where N 2 k-i(£) is as in (|2.7p but with u 2 k-2 replaced with u^k-i- Set e(t,x) = 
v(r(t), X(t)x) and y = \(t)x, r = r(t). Then, with A = ^, 

fte(t,r) = r'(t)(w r + AA" 1 ^) 

and 

(3.2) tt e(t, r) = r"(t)(9 r + X\~ l yd v )v + r'(t) 2 (d T + XX~ 1 yd y ) 2 v 
Set 

/•to 2 ! 

r(t) = / A(s rfs + -to v = -t~ v 
Jt v v 

so that t' (t) = A(t), and r"(i) = A(r)A(r) (we are writing A(r) instead of A(t(r))). 
Then (|3.ip can be rewritten as 

[{d T + XX^ydyfv + XX- 1 (d T + XX~ 1 yd y )v ~ Av - hW A v]{T : y) 

= A- 2 (r)[7V 2fc _ 1 ( e ) + e 2fe _ 1 ](<(r),A- 1 y) 
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We remark that the linear Schrodinger operator H = —A — hW A on L 2 (R 3 ) has 
at least one negative eigenvalue as well as a zero energy eigenvalue and resonance. 
The negative spectrum renders the linear evolution in (|3 . 3[) exponentially unstable. 
To address this further, we switch to the radial variable R = Xr from the previous 
section. Thus, (13.31) is the same as 



[(d T + XX~ X Rd R ) 2 v + XX- 1 (d T + XX-'Rd^v - v RR - \v R - 5W 4 v}{t, R) 

= A- 2 (r)[^ 2fe _i( £ ) + e 2 fc-i](t(r), X~ 1 R) 
or, with the new dependent variable e(r, R) := Rv(r, R), 

(d T + XX~ 1 Rd R ) 2 e - AA _1 (d T + XX~ 1 Rd R )i + Ce 
(,) ' 4 ' = X- 2 R{N 2k -i(R- 1 e) + e 2k - 1 ] 

where 

C = -dun - 5W i (R) on L 2 (0,oo) 
with a Dirichlet boundary condition at R = 0. Let 

0(t) := AA- X (r) = V := d T + P(r)Rd R 

TV 

and rewrite (13.41) in the form 

(3.5) V 2 e - (3{t)V£ + Ce = f 

To solve this equation we need precise spectral information on the operator C. 

4. The spectral and scattering theory of the linearized operator 
Definition 4.1. Let ^ 

C:= ~ dRR ~0TR 2 JW 2 

be the half-line operator on L 2 (0,oo) with a Dirichlet condition at R = 0. It is 
self-adjoint on the domain 

Dom(£) = {/eL 2 ((0,co)) : /, /' e AC([0, R}) V R, /(0) = 0, /" 6 L 2 ((0,oo))} 

Note that C<fi = where 

cj)(R) := 2Rdx ^ \*W(XR) = R(l - R 2 /3)(l + R 2 /3)-^ 

This means that C has a resonance at zero energy. Since <j) has a single positive zero, 
it follows from oscillation theory, see [3J, that there is an unique simple negative 
eigenvalue which we denote by Thus, there is cj>d € L 2 (0,oo) n C°°([0, oo)), 
decaying exponentially, and with (j>d(R) > for R > but (^d(O) = so that 
£-4>& — £,d4>d- We also assume that H^dlla = 1- Clearly, C has no other eigenvalues 
or resonances. 

Lemma 4.2. The spectrum of C equals 

spec(£) = {£ d }U[0,oo). 

The positive spectrum is purely absolutely continuous, and ^ is a simple eigenvalue 
with eigenfunction that we denote by 4>d- Moreover, C has a resonance at zero. 
In fact, C(/)o = with (f>a(R) = R(l - R 2 /3)(l + R 2 /3)'i. Finally, C is in the 
limit-point case at infinity. 
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As usual, see Marchenko [12] or Section 2 of [4], we introduce the standard 
fundamental system of solutions (f>(R, z) and 9(R, z) for all z G C of Cy = zy with 
the boundary conditions 

0(0, z) = 9'(0, z) = 0, <f>'(0, z) = 9(0, z) = 1 

so that in particular 

W(6(;Z),<t>(;Z)) = l 

These functions are entire in z. Note that 4>(R, 0) = 4>o(R) from above. Further- 
more, 

9 (R) := 9(R,0) = (1 - 2R 2 + i? 4 /9)(l + i? 2 /3)"i 
The Weyl-Titchmarsh function m(z) is uniquely defined by 

(4.1) V+0) z ) : = z ) + ™>(z)4>(;z) £ L 2 (0, oo) V Imz > 

The solution ip + is referred to as the Weyl-Titchmarsh solution. Then one has the 

following, see H): 

Proposition 4.3. The function m can be analytically continued to C\spec(£) and 
it is a Herglotz function. For each R>0, ip+(R, z) and i/j + (R,z) are analytic on 
C\ spec(£). The spectral measure of C is 

dp = S Sd + p(0 = - Im m(C + t0) 

7T 

in the following sense: the distorted Fourier transform defined as 

^ : / — ► / 

poo pb 

/(&)=/ MR)f(R)dR, /(0=lim / <t>(R,0f(R)dR, £>0 
is a unitary operator from L 2 (K. + ) £o L 2 ({^} U R + , p) and its inverse is given by 

t- 1 : f — ► /(i?) = ffa)<t> d (R) + lim /" W 0/(0 p(£) d£ 

Here lim refers to the L 2 (R + ,p), respectively the L 2 (R + ), limit. 

In the sequel we view the Fourier transform as a vector-valued map 

7(60 



f 

1 /(•) 

The Weyl-Titchmarsh solutions are scalar multiples of the Jost solutions f+(R, z) 
which are determined via the condition that 

(4.2) Cf + (;z) = zf+(;z), f+(R,z)~e i ^ R BB R^>00 

where Imz > 0, Im^/z > 0. They are solutions to the integral equation, with 
V(R) = -5(l + R 2 /3)-\ 

f + (R, z) = + ^ Sin ^ {R '- R)) V{R')f + {R!, z) dR> 

Jr V z 

The functions ip+(R,z) have well-defined limits as z — > £ + iO when £ > 0. In 

i 

particular, tp + (R,£) ~ c (£)e l ^ R as R — » oo. The constant c (£) is determined 
from the Wronskian condition W(?/'+(-, £)> 0) = 1- Once we have determined 
c o(£) we find m(£ + «0) from the Wronskian relation 

m(£ + »O) = W(0(-,0^+(-.£ + iO)) 
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We now give an asymptotic expansion of our fundamental system for small z. 

Proposition 4.4. For any z€C the fundamental system <j)(R,z), 6{R,z) admits 
absolutely convergent asymptotic expansions 

oo 

<f>(R,z) = MR) + R- 1 J2( R2z ) j( t>AR 2 ) 

oo 

0(R, z) = O (R) + J2( R2z ) j0 j( R2 ) 

i=i 

where the functions cfrj, 9j are holomorphic in O = {ti 6 C : Reu > —5} and 
satisfy the bounds 

& 
(~)! 



(4-3) |0 3 -(u)| < —- — \u\(u)-*, 



(4.4) <__( u ) 5 , ue n 



Furthermore, 
(4.5) ct>x{u) 



-\u{l + o(l)) as 

(l + o(l)) as moo 

Proof. We begin with <fi. We formally write 



<f>(R,z) = R- l Y,^fA R l MR) = R4>o(R) 

3=0 

This becomes rigorous once we verify the convergence of the series in any reasonable 
sense. The functions fj should solve 

CiR- 1 ^) = R^fj-i, fj(Q) = /j(0) = 
where we have set /_i = 0. The forward fundamental solution for C is 

h(r,r') = (MR)W) - MR')HR))1[r>r>] 

Hence we have the iterative relation 

f R R r i 
= / 157 MR)W) ~ MR'Wo(R) fj-i(R')dR', f (R) = RMR) 
Jo ft 1 J 

Using the expressions for tfio, 9q we rewrite this as 

fj(R) = 

" R R 2 {\ - i? 2 /3)(l - 2R' 2 + i?' 4 /9) - RR'(1 - i?' 2 /3)(l - 2R 2 + i? 4 /9) . , D „ , D , 

5 -. fi—\{R ) dR 

R'{l + R 2 / , Z)i{l + R' 2 /Z)i J 

It is clear from this integral that each fj(R) is an analytic function of R 2 provided 
Rei? 2 > —1. Moreover, fj(R 2 ) vanishes like i? 2 ^ +1 ^ around R = 0, and grows 
like i? 2j+1 as R — > 00. The bounds in (14. 3p are a quantitative version thereof, and 
proved by induction. 

For 6(R, z) we make the ansatz 

00 

6(R, z) = J2 z J 9j(R), 9o(R) = o (R) 
3=0 
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where the functions gj solve 



£9j=9j-i, .9-1=0 



Thus, the iterative relation is 



MR)0o(R') - MR'Wo(R) gj-i(R') dfi, g (R) = (R) 



The analyticity is the same as for the fj(R), and gj(R) vanishes like R 2j around 
R = 0, whereas the growth is R 2j+1 as R — ■> oo, as claimed. 

Finally, the leading order of 4>\ (u) is found by solving for the coefficients c\ , C2 

in 

(-d RR -5(l + R 2 /3)- 2 ) Cl R 3 (l + o(l)) = R(l + o(l)) as R -► 
(-<9fli?,-5(l + i? 2 /3)- 2 )c 2 i? 2 (l + o(f)) =-V3(l + o(l)) as i? ^ oo, 
respectively. Thus, c\ = — g, c 2 = as claimed. □ 
Next, we turn to the asymptotic expansion of the Jost solutions. 

Proposition 4.5. For any £ > 0, f/ie Josi solution /+(■,£) as m (|4.2|) is of the 

form 

f + (R,Z) = e iR ^o-(R^,R), R 2 t>l 
where a admits the asymptotic series approximation 

oo 

a(q,R)^J2^^t( R ) 

3=0 

in the sense that for all integers jo > 0, and all indices a, f3, we have 

30 



(4.6) sup (R) 2 (Rd R ) a (qd q f a(q, R) - T q- j 4>f(R) 

R>0 1 



3=0 



< C, 



c*Mo q 



-i'o-l 



for all q > 1 . Here 

4'+ = l, 4t(R) 



as R — ► oo 
as fi^O 



zcii?" 2 + iO(R- 4 ) 
ic 2 R + iO{R 2 ) 

with some real constants c\ , c 2 . More generally, ipt (R) o,re smooth symbols of order 
—2 for j > 1, i.e., /or all k > 

sup(R) 2 \((R)d R ) k 4UR)\<oo 

R>0 

Finally, ipf(R) = 0{R : >) as 0. 
Proof. With the notation 

O -( 9)J R) = / + (i?,0e- i ^ 
we need to solve the conjugated equation 

5 



(4.7) 



-d RR - 2i^d R - 



(l + i? 2 /3) 2 

We look for a formal power series solving this equation 



a(R^,R) = 



(4.8) 



J=0 
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This yields a recurrence relation for the fj's, 
with /o = 1, which is solved by 

tm - - 2 d R f^{R) - - 2 J r jYTww fj - l{R,) dR ' 

Then fi(R) is smooth for all R > with /i(-R) = -5^ + iO(R- 5 ) as i? -> 00 and 
= ic + iO(R) around # = 0. More generally, = iW(R- j ~ 2 ) as i? -> 00 

and fj(R) — VO{X) as i? — > 0. Differentiating these functions leads to symbol-type 
behavior as R — » 00. Defining ip~j~(R) :— R 3 fj(R) yields the desired bounds. 

To finish the proof, we first construct an approximate sum, i.e., a function 
cr ap (q, R) with the property that for each jo > we have 

(4.9) \(Rd R ) a (qd q f[a ap (q,R) - f>-^t(i?)] | < c a)/3j0 (i?)- 2 ^' 0-1 

i=0 

The construction of (J a p{q, R) is standard in symbol calculus; we set 

00 

a ap (q,R) :=5^g-^t(Ji)x(?«i) 
i=o 

where #j — > sufficiently fast and x is a cut-off function which vanishes around zero 
and is equal to one for large arguments. The bound (|4.9[) implies that cr ap (i?^5, i?) 
is a good approximate solution for (|4.7|1 at infinity, namely the error 

e(i?^,i?) = (-fcut - 2i^d R - (1 + ^ 2/3)2 ) ct «p(^0 
satisfies for all indices a, (3, j 

\(Rd R r(qd q Y j e(q,R)\ < c aM {R)S~ j 
To conclude the proof it remains to solve the equation for the difference <7\ = 



-d RR - 2i^d R - (1 + fl2/3)2 j <Ji(R^,R) = e(R£*,R) 

with zero Cauchy data at infinity. We claim that the solution <j\ satisfies 

KRdnnqdgfa^R)] < c aM q-i{R)-\ j > 2 

Note that this finishes the proof by defining a = a ap — o\. A change of vari- 
able allows us to switch from the pair of operators (Rd R ,qd q ) to (Rd R ,£,d^) with 
comparable bounds. We rewrite the above equation as a first order system for 
{vi,v 2 ) = {(J 1 ,Rd R cj 1 ): 

., 1 Vx \ j R^ 1 \ f Vl \ ( 

Cr 



V2 J \ (l+R 2 /3) 2 

Then we have 



R- 1 - 2i£* j\v 2 )-\Re 

d 



-j-\ v f>-R-i\v\ 2 -R\v\\e 
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which gives 
and by Gronwall 



^vl^-CiR-'M+Rlel) 



POO / T}/ \ C 

\v(R)\< J r (^J R'\e(R>)\dR> 
Then for large j we have 

(4.10) |e| < C i R~ 3 {R)~ i M < C i Rr°{R)~ 2 < q~ j (R}~ 2 

To estimate derivatives of v we commute them with the operator. For derivatives 
with respect to R we have 

d R (Rd R ) ( M - ( U 5K , \ , ) (Rd R ) ( V ; 




I 

— 10J?.(l--R 2 /3) 

V (l + fi 2 /3) 3 

But the right-hand side is bounded by from the previous step and the hy- 

pothesis on e, therefore as above ROrv is bounded by Rr' J . 
We argue similarly for the £ derivatives. We have 

. u 2 



Z8 6 (Re) 

The only difference is in the first term on the right, for which we write £5 = 
and we use the decay property of v with j replaced by j + 1 : 

\^V2\ < ^g-'- 1 < R-\~\ m(Re)\ < R-V 




as desired. Finally, higher order derivatives are estimated by induction using the 
above arguments at each step. □ 

Next, we describe the spectral measure of C. Due to the resonance at zero 
energy, the spectral density becomes unbounded for small £. In what follows, 
f-{R,0 :=f + (R,-0 = U(R,0- 
Lemma 4.6. For all £ > there is a(£) ^ so that 



<j>(R,Z)=a(Z)f + (R,Z) + a(Of-(R,0 

with 

i few ~ J 1 as £ -> 

[ £ 2 a S £ — > oo 

wifft symbol type behavior of all its derivatives. The density p(£) of the spectral 
measure satisfies 

' £-2 as £^0 
£2 as £ — > 00 
with symbol type behavior of all derivatives. 



p(0 
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Proof. By inspection, 

a{0 = w(f + U)J-U)) = Zt& w M->Q>U->® 

By the preceding asymptotic analysis we can evaluate 
with some small fixed e > to conclude that 



£2 as £^0 
as £ — > oo 



with the corresponding upper bound on the derivatives. This yields the desired 
upper bound on |a(£)|. To obtain the lower bound, we proceed as follows. First, 
observe that 

lm(/ + (i?,0/i(#,0) = -^ 

Second, it follows that 

Im (f_(R, QW(f+(R, 0, 4>(R, 0)) = 

so that 

|Ty(/+(fl,0,W0)l>^ | l ^ ( ( ^ 

^From our asymptotic analysis, again at R — e£,~^ , 

W(R,0\ > / 1 as £^0 



~U" 5 as £^oo 

which leads to the claimed lower bound on |W(</>(-,£), 
The Weyl solution 

V>+ {R, £ + iO) = + m(f + iO)(j)(R, £) 
satisfies tp + (-,£ + iO) = co (£)/+(■, £). Solving for c (£) yields 

m ^ ] w(/ + (-,o,0(-,O) |iy(/ + (-,o,0(-,O)l 2 

Since 

/ + W) = -^,^(/ + (',9,9W)) + «WW/ + (',0>^6) 

implies that 

-2i£* = 0, /- (-, 0) = -2* Im [W(6(; £),/+(•, (-, 0, 0)] 

we conclude that 

(4.11) ,(0 - i ImmK + ffl) - - 

The denominator was estimated above, leading to the desired bound on the spectral 
density. □ 
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5. The transference identity 

Returning to the radiation part e in ()3.5p . the idea is to expand it in terms of 
the generalized Fourier basis <f>(R,£) from Proposition ^. 31 i.e., write 

poo 

i(r,R) = x (t)MR) + / x(t,OHR,Op(0^ 
Jo 

and deduce a transport equation for the Fourier coefficients (xq(t),x(t,£)). The 
main difficulty in doing this is caused by the operator R8r which is not diagonal in 
the Fourier basis. We re-express this derivative in terms of the derivative 2£c>£. We 
refer to this procedure, which involves a certain error operator, as the transference 
identity since it allows us to transfer derivatives from R to £. We define the error 
operator K, by 

(5.1) Rd^u = -2£d 6 u + JCu 

where / = Tf is the "distorted Fourier transform" from Proposition 14.31 and the 
operator —2^8^ acts only on the continuous part of the spectrum. Apriori we have 

AC : Cg°({&} U (0, oo)) -> C°°({&} U (0, oo)) 

Splitting functions on spec(£) into a discrete and continuous component we obtain 
a matrix representation for /C, 



K. 



K*cd fccc 



Using the expressions for the direct and inverse Fourier transform in Proposition ^. 31 
we obtain 



K-dd — 


[R8r 


MR), MR)) 2 

1 L R 


K-dcf = 


a: 


f(0Rd R <t>(R,0p(0<%, MR)) , 

' L ~R 


ICcd(v) = 


(r8r 


MHR,v)) r2 

' L R 




<f 


f(£)Rd R <KR,t)p{t)dt, <f>(R,r))) 2 

' L R 


■ ( 


(jT 


2^f(OHR,Op(OdC,^(R,v)) T2 

1 L R 



ICdd = -L JCdcf = - I f(0K d (0p(0d^ lC cd (v) = K d ( V ) 



(5.2) 

Integrating by parts with respect to R in the first two relations we obtain 
1 

2' 

where 

K d {r,) = {Rd R M4>{R,r } ) 
Integrating by parts with respect to £ in (|5.2[) yields 

JCccfiv) = ( J KOlRdn - 2£dt]<l>(R,t)p(t)dt,<l>(R,Ti) 

A + "m)^ 
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where the scalar product is to be interpreted in the principal value sense with 

/ecg°((o,oo)). 

In this section, we study the boundedness properties of the operator JC. We 
begin with a description of the function and of the kernel Kq(ji, £) of K, cc . 

Theorem 5.1. a) The operator IC CC can be written as 

(5.4) k; cc = -(^ + ^M)^- 7? )+/Co 

where the operator K-o has a kernel Ko(r], £) of the forn^ 
(5-5) Ko (r,,0 = £ ^F(tv) 

with a symmetric function F(£,rj) of class C 2 in (0, oo) x (0,oo) satisfying the 
bounds 

\d e F(£,ri)\ + \d v F(Z,r,)\ < 
sup Wd*F(t,T,)\< 

j+fc=2 

where N is an arbitrary large integer. 

b) The function is smooth and rapidly decaying at infinity. 

Proof. We first establish the off-diagonal behavior of /C cc , and later return to the 
issue of identifying the ^-measure that sits on the diagonal. We begin with (|5.3p 
with / e C£°((0,oo)). The integral 

p OO 

<R) = / f(0[Rd R - 2Zd s ]<P(R,0p(0 d£ 
Jo 

behaves like R at and is a Schwartz function at infinity. The second factor (j)(R, rf) 
in (|5.3[) also decays like R at but at infinity it is only bounded with bounded 
derivatives. Then the following integration by parts is justified: 

vlCccf(v) = (u,£4>(R,vj) = (Cu,ct>{R^)) L2 









+ ?/ 






-?? < 1 


(t- 




_1 (n 








-?? > 1 








i 




C- 


fry < 1 


(e- 


\-vY 


-ni- 




-r)i\ 


)" JV £- 


1-77 > 1 






ce + 








-T) < 1 


«- 


br})' 


" 2 (H 








-77 > 1 



Moreover, 



/» CO />CO 

£u=/ f(0[C,Rd R ]<l>(R,0p(0(%+ f(0(Rd R -2^d e )^(R,0p(0^ 
Jo Jo 

I /"CO 



with the commutator 



2 The kernel below is interpreted in the principal value sense 



26 



J. KRIEGER, W. SCHLAG, AND D. TATARU 



Thus, 

Jo 

Hence we obtain 

V^ccfiv) ~ £cc(£/)(f7) 

I ■ - double integral on the right-hand side is absolutely convergent, therefore we 
can change the order of integration to obtain 

(V - OK (V, = P (0(U(R)<P(R, £), <p(R, v) 



mO(R9R-2^)4>(R,Op(0 dt 
f(0U(R)4>(R,0p(0dC,cl>(R,v)) 



This leads to the representation in (15. 5ft when £ ^ 77 with 
F(£,r)) = (u(R)4,(R,0,<t>(R,r))] 



It remains to study its size and regularity. First, due to our pointwise bound from 
the previous section, 



(5.6) 



sup|0(iU) 

_R>0 

\Rd R <f>(R,0 

\d 6 <P(R,0 

\d^(R,0 

\d!<t>(R,0 
\dl<p(R,0 

we always have the estimates 



< 



< R V£ > 1 

< min(i?^- 1 ,i? 3 ) 

< mm(RCKR 2 ) 

< mm{R 2 cKR 5 ) 
<mm(R 2 C\R 4 ) 



V£ > 1 
VO < £ < 1 

v<e > 1 

V < £ < 1 



(5.7) 



\%F(Z,ri) 



<rV x 

< 4 2^2 

< ^ 2 f] 2 



V£ > 1, V > 1 

> 1, 77 > 1 

> 1, 77 > 1 



They are only useful when £ and 77 are very close. To improve on them, we consider 
two cases: 

Case 1: 1 < £ + 77. To capture the cancelations when £ and 77 are separated we 
resort to another integration by parts, 
(5.8) 

vFfotj) = (u(R)4>(R,0,jCct>(R,v)) = ([C,U(R)}<p(R,0,(b(R,v)) +£F(Z,v) 



Hence, evaluating the commutator, 

(5.9) (T)-tMt,v) =-((2U R d R + U RR )<t>(R, 0,HR,V) 



Since U R (0) = it follows that (2U R d R + U RR )<f>(R,£) has the same behavior as 
4>(R, £) in the first region. Then we can repeat the argument above to obtain 

(V - 2 F(^v) = -([C,2U R d R + U RR ]<f>(R,0,<l>(R,ri)) 
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The second commutator has the form, with V(R) := — 5(1 + i? 2 /3) 2 , 

[C, 2U R d R + U RR ] = 4[/ flfl £ - Wounds. ~ U RRRR - 2U R V R - 4tf. RH V r 
Since V(R), U(R) are even, this leads to 

(v - 2 Ht v) = ((c/ odd (i?)a fl + [/—(i?) + tu even (R))<j>(R, v) 

where by JJ odd , respectively JJ even , we have generically denoted odd, respectively 
even, nonsingular rational functions with good decay at infinity. Inductively, one 
now verifies the identity 

fc-l k 

in , (v-0 2k mv) = ((E^'^f dR + J2^ en )^,0,H;V)) 
(5-10) j=0 1=0 

(R)mf d (R)\ + \U e k T n (R)\ < {R)- 4 - 2k Vj,£ 
By means of the pointwise bounds on (f> and 8r4> from (|5.6p we infer from this that 

( t) k -Hv)-* 

(v - 2k 

Combining this estimate with (|5.7|) yields, for arbitrary N, 

£ (f + *7) _1 (l + l£* ~V*\)~ N provided f + Tj>l, 

as claimed. For the derivatives of F we follow a similar procedure. If £ and ?y are 
comparable, then from (|5 .T[) . \d v F(£, rj)\ < Otherwise we differentiate with 

respect to 77 in (|5.10[) . This yields 

fc-l fc 

j=0 £=0 

Using also the bound on F from above we obtain 



fe-1 



respectively 



iw^)i< ^ ,, 2fc ^«i<e 



and 

(v-0 2 

which again yield the desired bounds. Finally, we consider the second order deriva- 
tives with respect to £ and r\. For £ and r\ close we again use the bound from (|5.7p . 
Otherwise we differentiate twice in (|5.10[) and continue as before. We note that it is 
important here that the decay of U£ dd and U^ en improves with k. This is because 
the second order derivative bound at has a sizeable growth at infinity which has 
to be canceled, 

\%<KR,0)\ «i? 4 
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Case 2: £, 77 <C 1. First, we note that -F(0, 0) = 0. This can be verified by direct 
integration, and is heuristically justified by the fact that U = [C, ROr) — 2£. The 
pointwise bound 

sup \d e F^, v )\<i 

t,v>o 

follows by differentiating ()5.8|) and from the bound \d^<p(R,^)\ < R 2 , see (|5.6[) . 
To bound the second order derivatives of F we recall the pointwise bounds, for 

o<e<i, 



\dU(R,0\ 



< 



i? 2j R < 



3 =0,1,2 



If rj < £ < 1, then these bounds imply that 



IVK,l)|< (R)- i R i dR+ i (R)- i R 3 C^dR+ ( (iZ)" 2 ^)-* dR 

JO J£ 7 Jr; _ 2 

<[i+iog(e/»?)]r' 

The logarithm in the middle integral is an artefact and can be removed using the 
oscillatory nature of dt<f)(R, £) in the regime R 2 £ > 1 as provided by Proposition ^. 51 
and Lemma 14.61 Loosely speaking, this means integrating by parts using that 

fy<t>(R,0 ~ R^dRe 1 ^ for R 2 £_ > 1 and small f. Thus, actually 

\d iv m,v)\<r^ 

A similar computation yields, for £ < 1, 

_ i 

\d!F(Z,v)\< / (R)- 4 R i dR+ / ^Rr^R^dRZC* 

JO J{~2 

This bound is too weak when £ -C 77 < 1. In that case, we differentiate (|5.9|) to 
obtain 



(77 - OdlF(t 77) = 2%F(£, 77) + (dt<KR, 0, (2^ fl a fl + C/ flJi )0(i?, 77) 



which in turn yields 
(5.11) 

-)2 



(v-0%F(Z,r,) = J [2%F(C,C) 
Using also the bound 

\RdR^(R,C)\ <mm(R£-i,R 2 ) 
we can evaluate the inner product in (|5.1ip as follows: 



0,{2U R d R + U RR )d ( <l>(;Q) d( 



<)7 



,Z),(2U R d R + U RR )d C <l>(;{) 



ft 5 ri ' r°o 

< / (R)- 6 R i R 2 dR+ / i (R)- 6 R 4 R(-z dR+ / i (R)- 6 R 2 ^- 1 RC" dR 

JO JQ~^ 

<[i+iog(c/o]r' 



cr 3 
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< 77* 



The logarithm appearing in the middle integral can be removed as before exploiting 
cancellations. Thus, (I5.11|) is controlled by 

\(r,-OdlF(£,r,)\<\J\- h d{ 

Since 77 3> £ this yields 

\%F(£,r,)\ < V -i 

which concludes the analysis of the off-diagonal part of the kernel. 

Next, we extract the 5 measure that sits on the diagonal of the kernel of JC CC 
from the representation formula (|5.3p . see also (|5.4p . To do so, we can restrict £,77 
to a compact subset of (0,oo). This is convenient, as we then have the following 
asymptotics of <j>(R,£) for R£z ^> 1: 

cf>(R,0=2Re [a{0e m ^]+O(R- 2 ) 

(Rd R - 2^)cf>(R,0 = -4 Re [& (Qe* 1 ** ] + 0(R- 2 ) 

where the O(-) terms depend on the choice of the compact subset. The R~ 2 terms 
are integrable so they contribute a bounded kernel to the inner product in (|5.3[) . 
The same applies to the contribution of a bounded R region. Using the above ex- 
pansions, we conclude that the <5-measure contribution of the inner product in ([57 
can only come from one of the following integrals: 



(5.12) -4/ / f(Ox(R)^[^'(0<v)e iR(r2+r, ' 2) ]p(0^dR 

Jo Jo 

r°° r°° , 1 i. 

(5.13) -4/ / f($ x (R)B*[Sa\$a(rj)e iR W^]p(QdedR 



where x is a smooth cutoff function which equals near R = and 1 near R = 00. 
In all of the above integrals we can argue as in the proof of the classical Fourier 
inversion formula to change the order of integration. Integrating by parts in the first 
integral (|5. 12[) reveals that it cannot contribute a ^-measure. On the other hand, 
(|5.13|) contributes both a Hilbert transform type kernel as well as a <5-measure to 
K. By inspection, the 6 contribution is 

Re[£a'(£)a(r7)e iR «*- 7) *)]p(£) dR 

) 

= -87r^Re[ea / (Oo(0]p(0*tt-»?) 

^ p(0 - 

where we used that piO 1 = 47r£5| a (£)| 2 in the final step, see (|4.11[) . Combining 



8(t - 77) 



this with the <5-measure in (|5.3p yields (|5.4|) . 
b) Arguing as in part (a) we have 

KM - ^ 

For F we use the representation in (|5. 10|) with £ replaced by £d and replaced 
by <fid- The conclusion easily follows from pointwise bounds on <fi(-,T)) and its 
derivatives. □ 
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Next we consider the L 2 mapping properties for /C. We introduce the weighted 
L 2 spaces L 2 ' a of functions on spec(£) with norm 

/>oo 

(5-14) := |/(^)| 2 + / \f(0\ 2 (® 2a p(0<% 

Jo 

Then we have 

Proposition 5.2. a) The operators Kq, K, map 

/C : L 2 p ' a -> L 2 p ' a+1 / 2 , /C : L 2 p a -> L 2 p ' a . 
b) In addition, we have the commutator bound 

: L 2 p a ^L 2 p a 

with acting only on the continuous spectrum. Both statements hold for all 

Proof. We commence with the /Co part, a) The first property is equivalent to 
showing that the kernel 

pHv)(v) a+1/2 K (v,mr a p-ko ■■ l 2 (r+)^l 2 (r+) 

With the notation of the previous theorem, the kernel on the left-hand side is 

We first separate the diagonal and off-diagonal behavior of Kq, considering several 
cases. 

Case 1: (£,7?) e Q := [0,4] x [0,4]. 

We cover the unit interval with dyadic subintervals Ij — [2 3 ' -1 , 2?' +1 ]. We cover 
the diagonal with the union of squares 

2 

j=-oo 

and divide the kernel Kq into 

1qK = IahqKq + Iq\aKq 

Case 1(a): Here we show that the diagonal part IadqKq of Kq maps L 2 to L 2 . 
By orthogonality it suffices to restrict ourselves to a single square Ij x Ij . We recall 
the Tl theorem for Calderon-Zygmund operators, see page 293 in [T^]: suppose the 
kernel K(r], £) on M 2 defines an operator T : S — ► S' and has the following pointwise 
properties with some 7 € (0, 1] and a constant Co: 

(i) 1^(77,^)1 < Col^ - ^l" 1 

(ii) \K( V ,0 - K(rf,S)\ <C \v~ VV\Z - vr 1 ^ for all |r? - rf\ < \£ - t?|/2 

(in) \K(r,, - K(n, o\ < coie - mz - vr 1 - 1 for a u ^ - e\ <\c- W2 

If in addition T has the restricted I? boundedness property, i.e., for all r > 
and £0,% € R, ||T(oM>)|| 2 < C* r^ and ||T*(w^''°)ll2 < C r^ where w r '«°(£) = 
— £o)A) w hh a fixed bump- function ui, then T and T* are L 2 (R) bounded 
with an operator norm that only depends on Co- 

Within the square Ij x Ij , Theorem 15.11 shows that the kernel of Kq satisfies 
these properties with 7=1, and is thus bounded on L 2 . 
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Case 1(b): Consider now the off-diagonal part \q\aKq. In this region, by 
Theorem 15.11 

lA-ofaoi <(^r* 

which is a Hilbert-Schmidt kernel on Q and thus L 2 bounded. 

Case 2: (£,??) £ Q c . We cover the diagonal with the union of squares 



B = |J li x I 3 

3=1 



and divide the kernel Kq into 



1qcK = IbhqcKo + 1qc\ B Ko 
Case 2a: Here we consider the estimate on B. As in case la) above, we use 
Calderon-Zygmund theory. Evidently, |-ftTo(?7>£)l S l£ — ^l -1 on B by Theorem l5.ll 
To check (ii) and (iii), we differentiate Kq, It will suffice to consider the case where 
the derivative falls on F(^,rj). We distinguish two cases: if — 772 | < 1, then 
|£ — t]\ < £2 which implies that 

if, on the other hand, |£3 — 77^ | > 1, then 

^1^4 V|£-£'|<K- 7/1/2 
|^-77||^-»ji| |£"??| 2 

which proves property (iii) on B with 7 = A and by symmetry also (ii). The 
restricted L 2 property follows form the cancelation in the kernel and the previous 
bounds on the kernel. Hence, Kq is L 2 bounded on B. 

Case 2b: Finally, in the exterior region Q c \ B we have the bound, with arbi- 
trarily large N, 

\K (v,O\<a+O- N ^ + v)- N 
which is L 2 bounded by Schur's lemma. 

This concludes the proof of the first mapping property in part (a). The second 
one follows in a straightforward manner since Kd is rapidly decaying at 00. 

b) A direct computation shows that the kernel Kq"" 1 of the commutator [£c\ , Kq] 
is given by 

K° m (v,0 = (r,d v + Zd s )K (r,,Z) + Ko(ri,0 = 0-F com (H,r,) 
interpreted in the principal value sense and with F com given by 

pcom^ v) = t^-Fit, 7?) + + TfaWt, V ) 

By Theorem 15.11 this satisfies the same pointwise off-diagonal bounds as F. Near 
the diagonal the bounds for F com and its derivatives are wors^l than those for F by 
a factor of (1 = . Then the proof of the L 2 commutator bound for Kq is similar 
to the argument in part (a). 

The remaining part of the commutator [/C,£<9,e] involves 



^The one derivative loss can be avoided by a more careful analysis, but this does not seem 
necessary here. 
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(i) The commutator of the diagonal part of /C cc with £e\. This is the multipli- 
cation operator by 

pb 

which is bounded since p has symbol like behavior both at and at oo. 

(ii) The operator S,d^lC c d which is given by the bounded rapidly decreasing func- 
tion £0^(0- 

(iii) The operator JC^d^ given by 

!C dc Pd e f= K d {P)Pd^f{P)dP = - f(Od 6 m d (0)dP 
Jo Jo 

which is also bounded due to the properties of K~d- □ 

6. The final system of equations 

We now rewrite the main equation (|3.5p in terms of the Fourier transform. With 
T as in Proposition ^. 3[ and with j3 = AA _1 , 



which gives 

f(d T + !3RdRj 2 = (d T + /3(-2^ + £)) T 

= (d T - 2j3PdA 2 J^+20}c(d T - 20Pd^+ (3 2 (K. 2 + 2[£%,/C])JF 
Recall that 

/>oo 

e(r, R) = x {T)<f> d + / x(t, P)<j>(R, P)p(0 dP 
Jo 

This leads to a transport type equation for the Fourier transform 

X(t) = (x (t),x(t,O) 

of £ by applying T to (|3.5[) . It is convenient to write it as a system for the two 
components: 

1 o (d T -2ppdA 2 + p ) x 

(6 1) 

= /3(J - 2K)(d T - 2(3Pd i )x - /3 2 (/C 2 - K. + 2[£% K])X 



+ \~ 2 J r R(N 2k -i(R~ 1 J r ~ 1 X) + e 2k -i) 

where it is understood that 

N 2k „ 1 {R- 1 T- 1 X) + e 2k - 1 = (N 2k _ 1 (R- 1 T- 1 X) + e 2k _ 1 )(t(T),\- 1 (T)R) 

Note that N 2k -\ and e 2k -\ are only defined on R < r, but for the Fourier transform 
we need to extend them to all R - this will be described in the next section, but for 
the moment just take an arbitrary compactly supported extension with the same 
regularity. 

We treat problem (|6.ip iteratively, as a small perturbation of the linear equation 
governed by the operator on the left-hand side. For this we need to solve the 
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following uncoupled system consisting of an elliptic equation and a transport 
equation: 



d 2 T + Co 



(d T - 2/3(r)^ s ) 2 + t] x(r, £) = 6(r, £), 



(6.2) 



We want to obtain solutions to (|6.1[) which decay as t — * oo. For the first equation 
above this is achieved by using the standard fundamental solution Hq which has 
kernel 



-IWk-oi 



This means that up to homogeneous solutions of the form e ^ d ^ T the unique 
bounded solution to the elliptic equation is 

1 



Xd{r) 



/ e -l^l 5 l T - CT l b d (a) da 
Jo 



The second equation is restricted to the range £ > 0. Hence it is hyperbolic in 
nature, which means we can solve it backward in time, i.e., with zero Cauchy data 
at t = oo. We denote by H the backward fundamental solution for the operator 



[d T -2(3( T )^) +£ 
and by H(r,a) its kernel, 

/oo 
H{r,a)f{a) da 

Combining the two components we obtain a fundamental solution for the system, 

n = diag(H Q ,H) 

Then we look for a solution X to (|6.1[) as a solution to the fixed point problem 
X = h((3(I - 2K)(d T - 2j3£d^X - (3 2 {JC 2 - K + 2[^,K])X 



(6.3) 



\- 2 TR{N 2k -i{R' 1 T- 1 X) + e 2k -i) 



Remark 6.1. One can also freely add Ce~^ d ' ,T to the first component x d of X. 
Thus the fixed point argument yields in effect a one parameter family of solutions 
X depending on the parameter C. 

The mapping properties of H are described in the following result, which was 
proven in [TP] , 

Proposition 6.2. For any a > there exists some (large) constant C = C(a) so 
that the operator H(t, a) satisfies the bounds 

( 'er\ c 



(6.4) 
(6.5) 

uniformly in a > r 



\\H{T,0)\\ , .: 

(d T -2P(r)^)H(T,o) 



< 
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This leads us to introduce the spaces L°°' N L 2 p ' a with norm 
\\f\\ LX , NLr :=supr N \\f(T)\\ L2 , a 

P T>1 P 

Then the above proposition immediately allows us to draw the following conclusions: 
Corollary 6.3. Given a > 0, let N be large enough. Then 

\\Hb\\ L ~,N-* Lr +^ + I (dr - 2P{T)Zd^Hb\ LaatN _^ a < C Q 1||6|| 



with a constant Co that depends on a but does not depend on N . 



The small factor N 1 is crucial here for our argument to work. For Hq we have a 
stronger straightforward counterpart of the above result: 

Lemma 6.4. The operator Ho satisfies the bounds 

||-ffoMlz,°°.- N + II^t-HoMIl 00 '™ — GivllMU 00 '^ 

We note that while the constant on the right cannot be small, we no longer lose 
powers of t compared to the bounds for H . Hence for fixed ./V we can choose To 
depending on N so that we gain the smallness in the range r > r . Combining the 
last two results we obtain 

Proposition 6.5. Given a > 0, let N be large enough. Then there exists tq 
depending on N so that for t > tq we have 

\\nb\\ L ^ N _ 2L2 , a+l/2 + (d T -2f3(T)^dAnb <Coh\b\\ LOO , NL 2,« 

with a constant Co that depends on a but does not depend on N . 

On the other hand, the nonlinear operator N2k-i from (|6.ip has the mapping 
properties stated in Proposition 16.71 below. We first relate the spaces L 2 p ,a to the 
Sobolev spaces in R 3 . 

Lemma 6.6. Let a > 0. Then 

||x|| L 2, Q x \\Rr 1 T~ 1 x\\ H z<*(m?) 

Proof. For integer k we have 

k k 

\\x\\ L *.* x w&f-^h* ~ E wr-'&f-^Wl^) 



j=0 j=0 



k 



3=0 



But 



R^CR = -d 2 R - ^d R - 5W\R) 



where the first two terms can be recognized as the radial part of the three-dimensional 
Laplacian. Hence we get 

k 

NIl;-* x £ II (- A - ^w\R)yR- 1 T- 1 x\\ L 2 m 

3=0 
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Since W is bounded together with all its derivatives, the conclusion of the lemma 
follows for integer a. 

For noninteger a we use interpolation. First we consider the map 

x i — * R~ x T~ x x 

and obtain the bound 

\\R~ 1 f~ 1x \\h^(«?) ^ IMIz,2.° 

To obtain the reverse bound we use the map 

u i— ► J-RS(u) 

where S(u) stands for the spherical average of a function u in R 3 . □ 
Proposition 6.7. Assume that N is large enough and g < a < j. Then the map 

x -» X- 2 TR{N 2 k-i{R~ X J r ~ X x)) 
is locally Lipschitz from L°°' N - 2 L 2 p ' a+1/2 to L°°' N L 2 ' a . 
Proof. Using the lemma, it remains to prove that the map 

v ^ \~ 2 N 2k -i(v) 
is locally Lipschitz from L°°' N - 2 H 2a+x to L°° N H 2a . We recall that 

N 2 k-i(v) = 5( u tk-i -Ua)v + lOui^-i v 2 + 10u^_ 1 v 3 + 5u 2k -i v 4 + v 5 , 

see the comment following (|3.1[) , The time decay is trivially obtained for all but 
the first term, for which we need an additional step, where we pull out a factor of 
u 2k -i — uq. Using the regularity of u 2k -i given in Theorcm l2.ll we obtain 



A* 

u 2k -i-u Q e ^2 IS 2 (R,Q k -i), A~ 2 ("Uq - u 2 k-i) e t- 2 S 2 (R- 2 , Qfc-i) 

This indicates that two units of decay in r are gained. On the level of Sobolev 
spaces we argue as follows: since we are working with inhomogeneous Sobolev 
spaces, we can localize the above estimate to unit cubes, as the £ 2 summability for 
N 2 k-i(v) with respect to unit cubes is inherited from any of the v factors. But 
in any unit cube Q the coefficients u 2 k-i have at most (1 — a)~^~~ singularities 
(where a = j ~ -p) so that we can bound them in Sobolev spaces 

A 2 A 2 

\\u2k-i\\m+2«(Q) < dist (Q Q y \\U2k-i - uo\\m+^(Q) < ^dist(Q,0) 
where we used that a < K. Then it suffices to establish the quintilincar estimate 

jj2a+l _ jj2a+l _ jj2a+l jj2a + l jj2a + l ^ jj2a 

which in three space dimensions holds for a > i (a standard application of the 
fractional Leibnitz rule and Sobolev imbedding, see [22], page 105). □ 
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7. Conclusion 

We now prove Theorem 11.11 We first construct a blow-up solution inside the 
cone as follows. We begin with the approximate solution v,2k~i and the error e2k-i 
given by Theorem l2.1l inside the cone. We extend them outside the cone to functions 
having the same size and regularity, supported in r < 2t. Then the relation 

2 

&2k-i = (~dt + d 2 + ~d r )u2k-\ + ulk-i 
r 

is valid only inside the cone. 

Using Propositions 16.51 16.71 we iteratively find a solution 

X e L°°> N - 2 L 2 p > a+1 l 2 

for the equation (|6.3p for t < to sufficiently small and 

1 v 

Then we set 

v = R-^-'X G L°°< N - 2 H 2a+1 

and 

u = U 2 k-1 + v 
Given the derivation of (|6.3p . the function u solves 

2 2 
{-d 2 + d 2 r + -d r )u + u b = {-d 2 t +d 2 + -d r )u 2k -i + ~ e 2k -i 
r r 

which implies that the function u solves the nonlinear wave equation 

(-d? + a 2 + -d r )u + u 5 = o 

r 

inside the cone. 

The second part of the argument is to extend the above solution u to the exterior 
of the cone K = {0 < t < t , 0< |x| <£}so that the blow-up occurs only at the 
tip of the cone. For this we first observe that the above function u is close to u 
inside the cone and close to outside, namely 

lim / \V(u(t) - u (t))\ 2 + \u(t) -u (t)\ 6 dx = 

and 

limy \Vu(t)\ 2 + \u{t)\ 6 dx = 

Hence given 6 > we can choose to so that the two quantities above are less than <5 6 . 
We let w be the solution to the equation 

{-d 2 + d 2 + -d r )w + w 5 = 
r 

with initial data 

w(t ) = u(t ), w t (t ) = u t (t Q ) 

Due to the finite speed of propagation we conclude that w = u inside the cone. To 
conclude the proof of the theorem we will show that w cannot blow up outside the 
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cone before time 0. For this it suffices to prove that the energy outside the cone 
stays small, 

(7.1) / \Vw(t)\ 2 + \w(t)\ 6 dx < S 

see Q2], Q1- 

This is proved using energy conservation. The energy £(w(t)) is conserved in 
time. At time to we have 

£(w(t )) = £{u(t ))=S(W) + O(S) 

hence at time t we must have a similar relation. But the energy inside the cone is 
already close to this, so we obtain 

On the other hand, we have the Sobolev inequality 



< S 



M 6 dx < 



AT 




with a universal constant (independent of t) . Combining the two inequalities above 
we see that for each t there are two possibilities. Either we have 



or 



f \{ut + \Vu\Ydx<5 
I \{ul + \Vu\fdx>\ 



'at 

The first alternative holds at t = to. Then a continuity argument shows that it 
must hold at all t, since the above integral is a continuous function of t for as long 
as it stays small. 
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